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A new formula to calculate the transport coefficients of the causal dissipative hydrodynamics 
is derived by using the projection operator method (Mori-Zwanzig formalism) in [f|. This is an 
extension of the Green-Kubo-Nakano (GKN) formula to the case of non-Newtonian fluids, which 
is the essential factor to preserve the relativistic causality in relativistic dissipative hydrodynamics. 
This formula is the generalization of the GKN formula in the sense that it can reproduce the GKN 
formula in a certain limit. In this work, we extend the previous work so as to apply to more general 
situations. 
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I. INTRODUCTION 



Hydrodynamic models have been extensively applied to analyze the collective aspects of relativistic heavy-ion 
collisions. These analyses have mainly been done so far for ideal fluids [2J. The effect of dissipation (viscosity and 
heat conduction) to this problem has started only recently and it is less well understood yet. One of the reasons 
for this is that, besides technical difficulties in numerical implementations, a relativistic extension of the dissipative 
hydrodynamics is not trivial at all conceptually [H 0, H, la 0, B S EH EH- A naive covariant extension of the 
Navier-Stokes equation leads to the problem of relativistic acausality and instabilities of the theory [T3, EH, EH, E3| • 

An essential factor to solve this problem is to introduce a memory effect with a finite relaxation time in the definition 
of irreversible currents pi. ITU ITU fl(| . An important point here is that, with the presence of memory effects, the fluid 
becomes non-Newtonian, that is, the irreversible current is not simply proportional to the thermodynamical forces. 

This raises several serious questions in applying the causal dissipative hydrodynamics to various phenomena at 
relativistic energies. The crucial one, we will focus in this paper, is that we cannot use the Geen-Kubo-Nakano (GKN) 
formula to calculate the transport coefficients, because the derivation depends on the Newtonian property of the fluid. 
See the discussion in Appendix [XJ To obtain the transport coefficients of the causal dissipative hydrodynamics, a new 
formulation should be developed. 

One possible approach to obtain the transport coefficients in the presence of memory effects is the so-called projection 
operator method (POM). The POM was originally proposed to obtain master equations and generalized Langevin 
equations from microscopi c dy namics by implementin g sy stematic coarse-grainings in terms of projection operators 
for macroscopic variables [l8l . EH, [2(| [U [22|, [H, H3, Hal It i s a l so known that the POM is useful to obtain the 
microscopic expressions of various transport coefficients [l], HU, HH, H3|. In the POM, the transport coefficients are 
related to the memory function of the generalized Langevin equation. Except for trivial cases, it is very difficult to 
evaluate the memory function exactly and some appropriate approximations are needed. A common method is to 
neglect a part of the projection operator in the memory function (as explained later). In the following, we call such 
an approximation as the Q approximation. It is known that the formula for the transport coefficients obtained with 
the Q approximation in the POM are equivalent to those of the GKN formula pi H2, HI, H3|- It is further known 
that the coarse-grained equation of a conserved quantity obtained by the POM with the Q approximation becomes a 
usual diffusion equation HH H3, HH • That is, the use of the Q approximation in the POM leads to the behaviors 
of Newtonian fluids. 

Recently, one of the present authors discussed the coarse-graining procedure in the POM without using the Q 
approximation pi ITtI 12a] . There, it was shown that the equation for a conserved number becomes a telegraph-type 
equation when the Q approximation is not introduced 25| . Note that the telegraphic equation is derived when a 
memory effect is introduced in a diffusion equation [25j]. This indicates that we can apply this method to construct 
the causal dissipative hydrodynamics in the POM, defining the microscopic expressions of the transport coefficients 
for non-Newtonian fluids in a consistent manner. Following this idea, new formulas of the transport coefficients for 
the causal dissipative hydrodynamics have been derived This new formula differs in a essential way from those 
obtained using the GKN formula with the Newtonian case, although it can reproduce the GKN formula under a limit 
where the Q approximation is valid. 

In this paper, we present a more detailed version of the work of [f| and derive more general expressions of transport 
coefficients, in particular, the shear viscosity of the causal dissipative hydrodynamics. This paper is organized as 
follows. In section [XH for the sake of later convenience, we review briefly the projection operator method to derive 
the generalized Langevin equation. In section [TTT1 the so-called Mori projection operator is introduced. We calculate 
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explicitly the memory function in section IIVI This result is the generalization of the formula obtained in [l| and one 
of the main results of this paper. By using this general expression, we define the causal shear viscosity coefficient and 
the relaxation time in section |V] The result of this section is completely same as that of The relation between our 
formula and the GKN formula is discussed in section PVTl In section IVTTI we apply the the result to an exactly solvable 
model to confirm the validity of our exact expression of the memory function. In section [Villi we reinvestigate the 
result of the section IVl and propose another possible definition of the causal shear viscosity coefficient. The section 
IIXI is devoted to concluding remarks. 



II. PROJECTION OPERATOR METHOD 

It should be emphasized that the projection operator method (POM) was firstly proposed by Nakajima [HI , although 
it is often refereed as the Mori-Zwanzig formalism due to the extensive use and developments done by these authors. 
This approach has been studied so far in various contexts of physics and chemistry [2^, [23|, [24[ ■ In particular, Mori 
introduced the so-called Mori projection operator to describe the dynamics near thermal equilibrium and derived a 
generalized Langevin equation from microscopic models [20| . The generalized Langevin equation (the Mori equation) 
gives the basis of the various development of statistical physics. Kawasaki, for example, developed the mode coupling 
theory which describes the dynamical critical phenomena by using the technique of the generalized Langevin equation 
(30| . The mode coupling theory is recently used to discuss the glass dynamics [3l|. It is considered that the POM is 
a promising method to establish a new coarse-grained dynamics like the dynamical density functional theory [32. I33I . 
The formulation of the projection operator method has been polished up by several authors [HI, [H, [H, l37l. l38l. KM 
Eol . l4lj . The derivation discussed here is following [3^| . 

In a quantum mechanical system, the time evolution of an operator is governed by the Heisenberg equation of 
motion, 

jO(t) = i[H,0(t)] 

= iLO{t) (1) 

— *0(t) = e^-^Oito), (2) 

where L is the Liouville operator and to is an initial time at which we prepare an initial state. In the following, we 
set to = 0. We consider here an isolated system so that the Hamiltonian is independent of time. Note that Eq. @ 
is also valid for classical cases provided that the commutator of the Liouville operator is interpreted as the Poisson 
bracket. 

In order to derive coarse-grained equations such as hydrodynamical equation of motion from a microscopic theory, 
we should construct a closed system of equations expressed only by those variables with macroscopic properties 
of the system. However, the Heisenberg equation of motion contains the information not only of gross variables 
associated with macroscopic (hydrodynamic) time scales, but also of microscopic variables. In the POM, the latter 
variables are projected out by introducing an appropriate projection operator P (to be specified later). We denote 
its complementary operator by Q(= 1 — P). They should satisfy, 

P 2 = P, (3) 
PQ = QP = Q. (4) 

Here, the projection operators are time-independent. To describe real non-equilibrium processes, in general, the 
projection operator should be time-dependent. However, for the purpose of the present paper, simple time- independent 
projection operators are suffice as the definition of the transport coefficients of the relativistic dissipative fluid. 

From Eq. (J2J), one can see that the time dependence of operators is determined by e lLt . This operator obeys the 
following differential equations, 

Multiplying the operator Q from the right, we have 

j t e lLt Q = e lLt PiLQ + e tLt QiLQ. (6) 

Equation © can be solved for e lLt Q, 

e iLt Q = Qe lLQt + f dre lLT PiLQe iLQ{t - T) . (7) 
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Substituting Eq. ([7]) into the last term in Eq. ([5]) and operating O(0) from the right, we obtain the so-called 
time-convolution (TC) equation, 

a rt 

—Oft) = e tLt PiLO(0) + / dTe lL ^- T ^PiLQe iLQT iLO(0) 
at J 

+Qe iLQt iLO(0). (8) 

The first term on the r.h.s. of the equation is called the streaming term and usually corresponds to collective modes 
such as plasma wave, spin wave and so on. The second term is the memory term that causes dissipation. The third 
term represents the noise term. The second term and third terms are related through the fluctuation-dissipation 
theorem of second kind, which will be discussed later. 

Discussion of this section has been done in the Heisenberg picture and the generalized Langevin equation ([5]) is 
derived. We can develop the similar discussion in the Schrodinger picture and obtain master equations. As we pointed 
out, to discuss more complex non-equilibrium processes, we have to change the basis of the projection with time. Then, 
the projection operator is explicitly time-dependent. As for the operation of the time-dependent projection operator, 
see [40l | and references therein. 

It is also possible to derive another form of the generalized Langevin equation, which is called the time- 
convolutionless (TCL) equation. There are some cases where we cannot implement the Markov approximation in 
the TC equation. The (f> 4 theory is one of the examples, and the Markov equation is derived from the TCL equation. 
See for details. 



III. MORI PROJECTION OPERATOR 



In the above derivation of the TC equation, we have not specified the projection operator P. As a matter of fact, 
there are many possible projection operators that extract the slowly varying components from dynamics. Suppose 
that the macroscopic dynamics can be described by the time evolutions of n-gross variables. Then we have to define 
the projection operator to project any time evolution onto the space spanned by these n-gross variables. For example, 
in the case of usual hydrodynamics, the time evolutions are described by the dynamics of the energy density, velocity 
field and number density. Then we have five variables which form a complete set for hydrodynamics (two scalar fields 
and one vector field). 

Strictly speaking, there is no general criterion to prepare a complete set of gross variables. It is, however, suggested 
that there are three candidates for gross variables [20l [30l ] : (i) order parameters (if there exists phase transitions) , 
(ii) density variables of conserved quantities and (iii) their products. Once we find out a complete set of the gross 
variables, any macroscopic variables should be approximately given by a linear combination of these gross variables 
(see below). We can use the Mori projection operator to implement this coarse-graining. 

Let us represents a set of gross variables by a n-dimensional vector, 

A T = (Ax,A 2 ,---,A n ). (9) 
Then, the Mori projection operator P is defined as 



po = J2 c * a *> ( 10 ) 



i=l, 

for an arbitrary operator O, where the coefficient Cj is given by 



Cl = J2(0,A])-(AA%\ (11) 



The inner product is Kubo's canonical correlation 



[ ^Tr[p e^Xe-^n (12) 
Jo P 

where p = e~^ H /Ti\e~" H \ with the temperature /3 _1 . The inverse of the canonical correlation is defined by 

^2(A,A)T j 1 .(A j) A k ) = 5 i , k . (13) 



As for the physical meaning of the Mori projection operator, see, for example, [171 [2 
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IV. THE EXACT EXPRESSION OF THE MEMORY FUNCTION 

Substituting the Mori projection operator into Eq. the TC equation is reexpressed as follows; 

^ A(t) = iAA(t) - J* drS(r)A(t - r) + £(t), (14) 

where 

iA = £(^,4)(AAt)-/, (15) 

S«(t) = -6{t)Y,{iLQe iLQt iLA h A\){A,A^}, (16) 

fe 

= Qe iLQt iLAi. (17) 

The TC equation using the Mori projection operator is called the Mori equation. 

The important information of dissipation is given by the memory function S(i). As a matter of fact, transport 
coefficients are defined by it. However, the calculation of the memory function is not simple because the expression 
of the memory function has the projection operator Q (see Eq. (JTSJ). As a matter of f act, the projection operator Q 
is approximately replaced with 1 to estimate the memory function in many textbooks [22 . 123 . l24j , 

B(t) » -Q{t)Y,{iLe iL HLA h A\){A,A\j. (18) 

As mentioned in the introduction, we call this procedure the Q approximation. 

Recently, Okada et al. calculated the memory function of the Ising model and found that the memory function can 
be expressed in terms of the combination of the usual time correlation functions [52]. Afterwords, Koide applied the 
same procedure to microscopic models and discussed the effect of the Q approximation [13, [H| . The coarse-grained 
dynamics of conserved quantities of the model are, usually, considered to be given by the diffusion equation. As a 
matter of fact, when we apply the Q approximation, we can derive the diffusion equation in the model. When we do 
not apply the Q approximation, however, the coarse-grained dynamics is given by the telegraph-type equation instead 
of the diffusion equation [2o| . More interestingly, if the model has a conserved quantity, we can derive the sum rule 
associated with the conserved law. It was shown that the telegraph-type equation derived with the memory effect is 
consistent with the sum rule, while the diffusion equation (Q approximation) breaks it [25L [27|. 

The same idea is used to define the transport coefficients of the causal dissipative hydrodynamics by using the 
simplest Moriprojection operator that is defined with only one gross variable [![. In this section, we extend the 
discussion of |l[ to more complex cases where the Mori projection operator is defined with n-gross variables. 

To calculate the coarse-grained time-evolution operator, we introduce the following operator, 

B(t) = 1 + VH) n / dt x --- dt n L p {t l )L p {t 2 ) ■ ■ ■ L p (t n ), (19) 

n=l J JO 

(20) 

L p {t) = e~ iLt PLe iLt . (21) 



with 



Then, the matrix 3(f) is rewritten as 

Eii(t) = -6(t)(iLe iLt B(t)QiLAi, A\){Aj, A^y 1 



-0(t) [(X(t))« - (X(0)X(t)) t3 

-J2(-^ n / dtl ■ ■ ■ / dt n (it(t n )±(t n -i - t n ) ■ ■ ■ X(ti - t 2 )X(t - h)) l3 

n=l Jo Jo 

00 ft l-tn-l 

- Y,(- l T / dt x --- dt„(X(0)X(t n )X(t n ^ - t n ) ■ ■ ■ X(ti - t 2 )X(t - t!))i 

Jo Jo 



(22) 
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See Appendix [B] for detailed discussion. By using the Laplace transform, the above expression can be rewritten as 
the following simple form, 

S L (a) = -X L (a) 1 + X(0)X L (s) 1 , (23) 

where the functions X i (s) and X L (s) are given by the Laplace transform of the following correlation functions, 

Xy(i) = J2( iL Mt),Al)(A,A%\ 

k 

Xy(*) = ^ iL ) 2 Mt),4)^A%\ (24) 
fe 

This is the exact expression of the memory function without using the Q approximation. One can see that if we 
set n = 1, the expression ()23[) reproduces the result of the previous result of The expression of the transport 
coefficients are derived by employing approximations to this memory function, as will will see in the next section. 

V. SHEAR VISCOSITY OF CAUSAL DISSIPATIVE HYDRODYNAMICS IN N=l FORM 

We apply the formula to define the shear viscosity coefficient of causal dissipative hydrodynamics. For simplicity, 
we consider a particular case of shear flow, where the fluid velocity points in the x direction and varies spatially in 
the y direction [H, Then, the energy-momentum tensor obeys the following equation of motion, 

where ir a P is the traceless part of the energy-momentum tensor, 

n kl = (sty - i<5%) T« (i, j, k, I = 1, 2, 3) (26) 

To define the shear viscosity coefficient, we introduce the Fourier transform of T 0x (y, t), and set O(0) = T 0x (k y , 0). 
And following the usual derivation of the shear viscosity coefficient of the GKN formula, we choose T 0x (k y , 0) as a 
unique gross variable [H, HJ]. In this case, the Mori projection operator is defined by 

PO = (O,T Qx (-k y ,0))(T Qx (k y ,0),T 0x (-k y ,0))-\ (27) 

where O is an arbitrary operator. We will discuss later a more involved case where we need two gross variables to 
define the Mori projection operator. 
Then, the TC equation ([5]) is given by 



^-T 0x (k y ,t) = iA(k y )T 0x (k y ,t) - J dT~(k y , t - r)T 0x (k yi r) + £(fc„, t), 



(28) 



where 

iA(k y ) = (iLT Ox (k y ,0) : T Ox (-k y: 0))(T Ox (k y ,0) : T Ox (-k y: 0)y\ (29) 

H(M) = I Z L (k y ,s)e st ds, (30) 
2m J Br 

Z(k y ,t) = Qe lLQt iLT 0x (k y ,t). (31) 
The Laplace transform of the memory function S i (fcj / , s) is given by 

Z L (k y ,s) = -X L (k y ,s) 1 -, (32) 

1 + X L {k y , s) 

where 

X{k y ,t) - (iLT 0x (k y , t), T 0x (—k y , 0))(T Ox (k y , 0), T 0x (—k y , 0)) _1 , 
Xa{ky,t) = ((tL) 2 T 0x (k v ,t),T ax (~k v ,0))(T 0x (k y ,0),T ax (-~k v ,0))- 1 . (33) 
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In this derivation, we used 

X(0) = (iLT 0x (k y ,0), T 0x (~k y , 0)) = 0. (34) 

So far, everything is exact formally. Now we carry out the coarse-grainings of the time scale to break the time- 
reversal symmetry. For this purpose, first of all, we separate the memory function into the two terms as follows 

^-T° X (k y ,t) = - f dTn 2 (k y ,t-T)T° X (k y ,T)~ [ dT<S>(k y ,t-T)T° X (k y ,T). (35) 



uz JO JQ 

Here, we dropped the noise term. The frequency function and the renormalized memory function are defined by 

n 2 (k y ,t) = i J ^-Im[3 L (k y ,-iu + e)]e- iut , (36) 
Hky,t) = y^Re[S L (fc„-^ + e)]e- Mt , (37) 

respectively. 

To introduce the coarse-graining in time, we have to know the temporal behavior of the two functions. The behavior 
of the two functions have been investigated for some special cases, such as the chiral order parameter in the Nambu- 
Jona-Lasinio model fl7l . exactly solvable model of many harmonic oscillators [l7| and the non-relativistic model with 
a conserved density [251 ]. For all these cases, the two functions exhibit common properties; the frequency function 
converges to a finite value and the renormalized memory function vanishes at late time. Inspired by these examples, 
we introduce an important assumption that these features for the temporal behavior of the two functions are valid 
in general. That is, the renormalized memory function relaxes rapidly and vanishes at large t, while the frequency 
function converges to a finite value after short time evolution. In principle, the validity of the assumption should be 
checked for more general examples by implementing numerical calculations. Once we accept the above assumption, 
we may introduce the following ansatzs for the memory functions incorporating these basic features essentially [l7l ]: 

n 2 (ky,t) — » D ky ky 2 , <$>{ky,t) — » — 6 (t) , (38) 



where 



D ky = — y lim n 2 (k y ,t), (39) 

ky t-OO 

1 f°° 

= / dt<£(ky,t). (40) 

Thy JQ 

The factor k 2 is introduced for the later convenience (see Ea. (l43l) b The above ansatzs are shown to be consistent with 
the final value theorem of the Laplace transformation [ll. 1 171 l25j. That is, when the renormalized memory function 
converges to zero at late time, its Laplace transform $ L (fe y , s) should satisfy [43| . 

lim $(ky,t) = lims$ L (fc a ,,s) = 0. (41) 

t — >oo s — >0 

Similarly, for the frequency function, 

D k k 2 y = hm n 2 {k y ,t) = lim s(n 2 ) L {k y ,s). (42) 

Using these expressions, we have the equation for the energy momentum tensor component, 

d /"* 1 

^Ox/i i\ _ n l2 / J ..XfU ^ rriOx/ 



T^(k y ,t) = -D ky k 2 dTwu x (k y ) r> x (k y ,t). (43) 

Ol Jq T ky 

Here, we expressed the x-component of the fluid velocity as 

U*(ky)=T *(ky)/w, (44) 

where w is an enthalpy density [l[. 
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On the other hand, as pointed out before, the time evolution of the energy-momentum tensor in causal dissipative 
hydrodynamics is given by a kind of the telegraph equation . In particular, in the special case discussed here ( the 
fluid velocity points in the x direction and varies spatially in the y direction), the linearized equation of the causal 
dissipative hydrodynamics is given by the following telegraph equation 0, |Tl[ , 

Q2 i a NN n2 

w T '" + ^ T '" + 1 —w u ' = '>- (45) 

This equation defines the causal shear viscosity coefficient 77^^ and corresponding relaxation time r. By comparing 
Eq. (|4"3")l with Eq. (j4"5|) , we obtain the expression for the causal shear viscosity coefficient and the respective relaxation 
time as 

7} NN = ]xm o 2wD kv T. (46) 

which are the results obtained in In this derivation, it is assumed that the projection operator is defined by only 
one gross variable. We will reconsider this derivation in section IVIIH 



VI. GREEN-KUBO-NAKANO FORMULA 



We obtained an expression of the causal shear viscosity coefficient for non-Newtonian fluids. On the other hand, it 
is well-known that the shear viscosity coefficient of Newtonian fluids is given by the GKN formula. In this section, 
we discuss the relation between our formula and the GKN formula 

As was mentioned before, it is known that, when we apply the Q approximation, the GKN formula is reproduced 
in the pro jection operator method [4ij . When we apply the Q approximation, the memory function is given by 



E L (k y ,s) w -X L (k yi s). (48) 



22, 23, 24 



On the other hand, when the correlation function X L (k y , s) is very small, the memory function (|32[) is then expanded 
as follows, 



E (ky,s) = —X (k y ,s) + X (ky, s)X (k y , s) . (49) 

That is, the correlation function X L {k y , s) represents the correction to the Q approximation. 

As a matter of fact, in the Q approximation, we can derive the relativistic Navier-Stokes equation and then our 
formula reproduces the GKN formula of the shear viscosity coefficient. First of all, the correlation function X(k y ,i) 
is rewritten as 

X L (k,s) = J™ dte~ st 9(t)^ J d 3 xd 3 x 1 e- lkx k 2 ( 7 r^(x,t),7r^(x 1 ,0))(T te (x 1 ,0),r Oa; (O,0))- 1 

-1 POO p poo 

= ~Wf3Ja dt J d3xd3x ^ St ~ lkXk2 J Q ^( 7 r Q ' 3 (x,rV Q/3 (x 1 ,0)) I . et (T te (x 1 ,0),r te (O,0))^ 1 , (50) 

where 

<7T^(X, t)n a0 (x U S))ret = ~ ^([^(x, t), TT af3 (xi , s)]) eq . (51) 

In this derivation, we used the relation [45| 



(tt/xi/, TTpcr) = f A MP A„ CT - A^A^p — -A^A pa ) , (52) 



where L„ is a scalar function and 

A^ = g^u - u^u„, (53) 
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with it M the four- velocity of the fluid in the Landau frame. The correlation function X L (k, —ioj + e) is real in the 
low momentum limit. Then, the frequency function vanishes and the equation of the energy- momentum tensor (|43[) 
is given by the linearized relativistic Navier-Stokes equation, 

9 T 0x _ v NS k 2 T 0x = o (54) 

at y 
where the Navier-Stokes shear viscosity coefficient is 

km — = - V ns \l 2 . (55) 

By using the expression of Tk, the Navier-Stokes shear viscosity coefficient is expressed by using the time correlation 
function as follows, 

V NS = I" dtO(t)± J d 3 x 1 (^(x,t),7r Q/3 (x 1 ,0))(T te (x 1 ,0),T te (O,0))- 1 . (56) 

Except for the normalization factor (T 0a: (xi, 0), T 0x (0, 0)) _1 , this expression is nothing but the GKN formula of the 
shear viscosity coefficient 5EI . That is, our new formula can reproduce the result of the GKN formula when the 
correlation function X L {k yi s) disappears in the low momentum limit. In this sense, our formula is the generalization 
of the GKN formula. 

As was pointed out, the vanishing X L (k y ,s) corresponds to the Q approximation. So far, because the exact 
expression of the memory function (|23l) was not known, we could not discuss whether the Q approximation is applicable 
in the low momentum limit or not. Now the validity of the Q approximation can be quantitatively estimated by 
calculating the correlation function X L {k J , 1 s). In fact, it is already known that there are examples where the Q 
approximation cannot be applicable [13, Hfl 

VII. MODE-COUPLING THEORY OF DENSITY FLUCTUATIONS 

So far, we discussed the simplest case where the system has only one gross variable. In this section, we will consider 
a more complex case where we need two gross variables to define the Mori projection operator. Such a situation will 
occur, for example, in a glass dynamics [3l| . Glass is a high density system and a particle is thickly surrounded by 
other particles. The energy and momentum of the particles are continuously exchanged by collisions. However, it is 
difficult for particles to move away from the initial position because the space around has already occupied by others. 
This is called jamming. In a glass dynamics, we usually choose the gross variables as the fluctuations of density of 
particles and the corresponding current. 

We consider a classical TV-particle system, where the Hamiltonian is given by [3l[ 



Then, the corresponding Liouville operator is given by 



In this case, the Fourier transform of the fluctuations of the particle number density pk (t) is given by 

SpK(t) = Y / e lk ' Mt) -(27r)V(k), (59) 

i 

where p = N/V and V is the volume of our system. In this system, the number of particle is a conserved quantity 
and the fluctuations of the density should satisfy the equation of continuity, 

<5p k (i) = i|k|jk(*). (60) 

Here, wc define the current, 

3k = - V(k • Pi )e ik " r S (61) 
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where k = k/|k|. Then, we choose the set of the gross variables as follows, 

a-('£). m 

By substituting into Eq. ([8]), we obtain the evolution equation of A. We further multiply from the right and take 
the thermal expectation value. Then the evolution equation of the correlation function is given by, 



8_ 

Of 



C(t) = »AC(t) - / dr3(T)C(t - r), (63) 
Jo 



where 

C(t) = (A(i)At) eg 

0'k(*)fy-k}eg (jk(t)j-k)eq 



(64) 



Because of the equation of continuity (|60p. these four correlation functions are not independent. It should be noted 
that, in the usual discussion of the glass dynamics, we do not consider the thermal equilibrium environment discussed 
here and the calculation of the memory function is more involved. 
The coefficient of the streaming term is given by 

( ik \ 

iA = ik Q , (65) 

where the static structure factor is S(k) = 4f(p_k(0)j0k(0)). By using Eq. (|23|) , we found that the upper components 
of the matrix of the memory function vanish, 

B L (s)=( ° ° Y (66) 



J 21 



is) SUs) 



Because of A*) = 0, the noise term disappears. 

We will concentrate on the element in the lower column of the matrix. Then, we can obtain the following two 
equations, 

£ F ( fc '*) + ^(k) i? ( fc '*) + j Q dT(^ 21 {r){ik)F{k,t-T)+S 22 {r)^ t F{k,t-T)\ = 0. (67) 

8 3 k 2 8 /"* / 8 8 2 \ 

W^ t) + ^^)8l F ^ t) + L ^^WW/(M-r) +H22 (^ f (M-r)) = 0. (68) 

Here, we introduce the following function, 

1 

N 



F(k,t) = -(Sp k (t)5p^) eq . (69) 



From the consistency of the two equations, one can find that H 2 i(r) should vanishes. As a matter of fact, this is 
shown by using the following exact relation, 

= &(-*),£<)* - m) ■ (A, At)] tf . (70) 

This relates the memory function and the noise term and is called the fluctuation-dissipation theorem of second kind 
(26|. From this relation, one can show that H 2 i(£) also should disappear when Bi2(t) vanishes. The correlation 
functions, then, should satisfy the following relation, 

3Us) oc (X&(*) - [X(0)ir L ( s )] 21 )(l + X&(*)) - (X*( a ) - [X(0)X L ( S )] 22 )X^( S ) = 0. (71) 

By using this relation, we can simplify the remaining memory term, 

„L (e , ^ 2 ( S )-[X(0)X^( S )] 22 



1 + X£ 2 ( S ) 

s 2 F L (k, s) - sS(k) + k 2 F L (k, s)/(mPS(k)) 
sF L (k,s) -5(k) ' 



(72) 
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This is the result obtained by using our expression of the memory function (|23|) . 

In the case discussed here, however, all correlation functions are expressed by the unique correlation function F(k, t) 
and we do not need to use Eq. ([23]) to calculate the memory function. From Eq. (|67|) . the Laplace transform of the 
equation is 

s 2 F L (k,.s) -sF(k,0) -F(k,0) = - m ^ k) F L (k, S )-S^(s)(sF L (k,s)-F(k,0)). (73) 

The expression of the memory function that is obtained by solving the equation above is same as Eq. (|72p . This 
means the consistency of our formula. 



VIII. SHEAR VISCOSITY OF CAUSAL DISSIPATIVE HYDRODYNAMICS IN N=2 FORM 

In the previous derivation of the causal shear viscosity, we assumed that the macroscopic motion can be projected 
onto the space spanned by the unique gross variable T 0x (k y ). If this assumption is correct, the memory function 
converges to a constant rapidly and we can define the causal shear viscosity coefficient and the relaxation time as was 
done in section [Vj The memory functions are calculated so far for several examples and the behaviors are consistent 
with this assumption. 

However, there is another suggestion for the definition of the projection operator. In the wake of the discussion of 
the extended thermodynamics [7j, Ichiyanagi proposed that the Mori projection operator should be defined by using 
not only usual hydrodynamic variables but also the corresponding currents [46| , although any calculable formula was 
not given. In this section, we rederive the formula for the causal shear viscosity coefficient following his idea. 

The set of the gross variables are given by 



T° X (ky) 

T yx (k y ) 



(74) 



By substituting into Eq. we have 

8 



d T 0x (k Vl t) = -ik y Ty x (k y ,t), (75) 



§tT*°(ky,t) = -ik y R k T 0x (k y ,t) - £ drS 22 (T)Ty x (k v ,t-T), 



(76) 



where 



Rk y = {Ty%ky) 1 Ty X {~ky)){T^{ky)^{-ky))-\ (77) 

The memory function is given by 



X 2 L 2 (s) + zk y R k X^(s) 

" 22(S) " ~ i + x-( s) 
i + x£ 2 ( s ) 

with the Laplace transforms of the following correlations, 



(78) 



X(t) = 



(T 0x (k y ,t),T 0x (-k y )) (T°*(k y ,t),Ty*(-k y ))\ ( {T°*{k y ),T°*{-k y ))- x 

(Ty x {k y ,t),T 0x (-k y )) (Ty x (k y ,t),Ty x {-k y )) J \ o (T^fo),!^-^))- 1 



X(t) = ^X(t), (79) 

X(t) = |^X(t). (80) 

Here, we omitted the noise term. 

It should be noted that we still define Kubo's canonical correlation by Eq. (|12p . where the expectation value 
is calculated by the usual thermal equilibrium state. This is different from the original idea of Ichiyanagi and the 
extended thermodynamics, where the concept of the thermodynamic variables are extended and hence we have to use 
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a non-equilibrium state to calculate the expectation value. However, when we restrict ourselves to the non equilibrium 
states whose deviation from equilibrium is still small, then the expectation values can be evaluated at the equilibrium, 
since the effect of the non-equilibrium expectation should be higher order and we assume that it is negligible. 

Equation (|75[) is the equation of continuity. If we can derive the causal dissipative hydrodynamics from the 
Heisenberg equation of motion, Eq. (|76|) should be reduced to the telegraph equation, 



g-T yx (k y ,t) = _^-(^K(^,i)--T^(fc,,t). (81) 



It should be noted that when we combine this equation with the equation of continuity, we can reproduce Eq. (I45|) . 
To obtain the telegraph equation from Eq. (I76p . we assume that the memory function H22 is Markovian; the memory 
function quickly vanishes with time, 

B a2 (t) = —*(*), (82) 



where 



1 r 00 

— = / drS 22 . (83) 

Thy JO 

By substituting them into Eq. (|76[) . we obtain 

^-Ty*(k v ,t) = ~ik v wR ky u x {k y ,t) - —Ty x (k y ,t). (84) 

By comparison with the telegraph equation, we identify the causal shear viscosity coefficient and the relaxation 
time as follows, 



NN 



= lim 2R ky r. (86) 

ky^O V 



Note that these formula look different form those in the previous section, formula (|47[) and (|46[) . We know that the 
two approaches (n = 1 and n = 2) shown in this paper are completely equivalent if no approximation is introduced, 
that is, the coupled equation (|75|) and (|76|) gives exactly same result as Eq. (|28|) . We even checked the consistency of 
the two approaches by solving the coupled harmonic oscillator model which is exactly solvable [l7T | . 

To use the formula and (156")) , we have to calculate three correlation functions, while we need two correlation 
functions in the formula (I47|) and (|46|) . Thus we should usually use the formula (|47|) and (|46|) to estimate the causal 
shear viscosity coefficient in causal dissipative hydrodynamics. However, when the calculated memory function does 
not satisfy the condition (f3"8")> . we have to use the formula and (|56"j) . 



IX. CONCLUDING REMARKS 



In this paper, we derived the general expression of the memory function extending the result of |l[. By using the 
expression, we define the shear viscosity coefficient and the corresponding relaxation time of the causal dissipative 
hydrodynamics. Our formula is the generalization of the GKN formula because, when the Q approximation is justified 
in the low momentum limit, the GKN formula is reproduced. 

Phenomenologically, the causal hydrodynamics is derived by introducing the memory effect to the relation between 
irreversible currents and thermodynamic forces. Thus the vanishing relaxation time limit r — > corresponds to the 
limit of the Newtonian fluid and hence the causal dissipative hydrodynamics is reduced to the relativistic Navier-Stokes 
equation (the Landau-Lifshitz theory). Thus it is sometimes expected that the causal shear viscosity coefficient is 
still approximately given by the calculation of the GKN formula, when the relaxation time is not large. However, 
this expectation is not trivial. As was discussed in this paper, the new formula reduced to the GKN formula in the 
Q approximation. In this limit, as is shown in Eq. (|55| . the causal shear viscosity coefficient t] NN vanishes and the 
expression of the relaxation time r is reduced to that of the shear viscosity coefficient in the GKN formula. That is, 
what is approximately given by the GKN formula is not the causal shear viscosity coefficient but the relaxation time. 

By using the idea of the AdS/CFT (anti-de Sitter/conformal field theory) correspondence in the string theory, we 
can calculate the correlation function of the energy-momentum tensor in J\f — 4 supersymmetric Yang-Mills theory 
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[47| . From this result, we obtain rj/s — 1/(4tt) with s being the entropy density and many people expect that this gives 
the minimum of the shear viscosity coefficient of the relativistic fluid of quarks and gluons. It should be, however, 
noted that the r] here is not rj NN but 77 , that is, this discussion is true only for Newtonian fluids, because, to derive 
the result, the expression of the GKN formula of the shear viscosity coefficient is used. Thus, when we discuss the 
causal dissipative hydrodynamics, we cannot use this value as the limit of the causal shear viscosity coefficient. The 
lower bound of the shear viscosity coefficient may exist even for the causal hydrodynamics. This will be predicted by 
using our new formula instead of the GKN formula. 

There are several approaches to derive the relativistic hydrodynamics consistent with causality. However, as far 
as we know, the telegraph equation plays an essential role to solve the problem of acausality in all theories, and the 
difference of the theories comes from the non-linear terms. Thus the formula discussed here is applicable even for other 
causal dissipative hydrodynmaics, the Israel- Stewart theory See [ll[, for more discussions about the relationship 
between different theories The effect of non-linearity, in general, can change the coefficients of the linear terms. To 
discuss the effect of non-linearity to the transport coefficients, we have to consider the non-linear response (49[. In the 
projection operator method, this is implemented by generalizing the projection operator including non-linear terms. 
However, the quantitative effect has not been known so far. 

On the other hand, the telegraph equation may be not unique solution of the problem of acausality in hydrodynamcis. 
For example, there are different approaches to solve this problem in diffusion processes 48]. However, to our best 
knowledge, there is no formulation of causal dissipative hydrodynamics in these alternative scenarios. It should also 
worth mentioning that we have not so far encountered any problem in implementing numerical simulations of the 
causal dissipative hydrodynamics [tj [ll[ . 

It should be mentioned that the projection operator approach discussed here and the usual linear response theory 
do not have explicit Lorentz covariance, because we introduce a thermal equilibrium background. That is, when the 
transport coefficients of relativistic fluids are calculated, we assume the existence of a local rest frame where the 
dynamics of macroscopic quantities are determined in a non-relativistic way, together with an appropriate boundary 
condition. 



T. Koide acknowledges helpful discussions with D. Hou. This work is supported by CNPq and FAPERJ. 



APPENDIX A: THE GREEN-KUBO-NAKANO FORMULA 



In the appendix, we gives the short review of the GKN formula. As for the calculations of the GKN formula for 
relativistic fluid, see, for example, [50| and references therein. 

We consider the system whose Hamiltonian is given by if. By applying an external force, the total Hamiltonian is 
changed from H to H + H ex (t), with 

H ex {t) = —AF(t), (Al) 

where A is an operator and F(t) is the c-number external force. 

We consider the current J induced by the external force. From the linear response theory, we obtain 

(J) = I dsV(t- s)F(s), (A2) 

J — OG 

where the response function is given by 

*(t) = / d\(A(-iX)J(t)) eq . (A3) 
Jo 

This is the exact result in the sense of the linear approximation. This formula, also, is called the GKN formula. 
However, in particular, when we define transport coefficients of hydrodynamics, we do not use this expression. 

In these cases, first, we assume the linear relation between currents and the external force, J(t) = DgknF{£) with 
the transport coefficient Dqkn- The formula to define the expression Dqkn is the GKN formula which is discussed 
in this paper. For this, we see that we should ignore the memory effect (time-convolution integral) in Eq. (|A2[) . 

/>OC 

(J) pa / ds<S>(s)F{t) (A4) 
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Then the GKN formula is 

D GKN = / dstf(s). (A5) 
Jo 

Thus this formula is applicable only when there is a proportional relation between a current and a force, like 
Newtonian fluids. This is the reason why we cannot use the GKN formula to calculate the transport coefficients of 
the causal dissipative hydrodynamics. 

In principle, it is possible to derive the transport coefficients of the causal dissipative hydrodynamics from Eq. 
(|A2|) . Instead of J(t) — DqknF^), we assume the following telegraph equation, 

d t J(t) = - — J(t) + —F(t). (A6) 

TR T R 

From Eq. (|A2|> . we can derive the following equation, 

poo 

d t J(t) = V(0)F(t) + / dsd a V(s)F(t). (A7) 
Jo 

In the second term, we ignore the time-convolution integral. We further assume the GKN formula to reexpress the 
first term. Then we finally obtain 

*(0) f°° 

d t J(t) = jr^J(t) + / dsd a *(a)F(t). (A8) 

L>GKN Jo 

By comparing this equation with Eq. (|A6|) . we can derive the expression of D and tr. 

Exactly speaking, we considered here the current induced by the external force. However, the shear viscosity 
is induced not by the external force but by the difference of the boundary conditions. Thus the discussion is not 
applicable to the problems discussed in this paper. 

APPENDIX B: THE DERIVATION OF EQ. (22) 

In this appendix, we derive Eq. (|2"2"j). By using Eqs. I|16p and (fTT))) . we obtain 

Ey(t) - -e^Y^iiLe^B^QiLA^Al^A,^)^. (Bl) 

k 

The first three terms can be calculated as follows, 

{iLe iLt QiLAi, A])(Aj, A])' 1 = X y (t) - [X(0)X(t)]^ (B2) 

(iLe iLt (-i) ( dse~ iLs PLe iLs QiLA u A^)(Aj,A^)~ 1 = - / ds[[X(s)X(t - s)]^ 
Jo Jo 

-[X(0)X( s )X(f- s )] y ], (B3) 



{iLe iLt f ds ie - iLsi PiLe iLsi f ' ds 2 e" lLs2 PiLe lLs2 QiLA l ,A])(A j ,A])- 1 
Jo Jo 

= I dsx I ds 2 [[X(s 2 )(X(si - s 2 )X(i - si)]y 
Jo Jo 

-[X(0)X(s 2 )X( Sl - s 2 )X(i - (B4) 



In short, the n-th order term is given by 



(iLe iU (-i) n ( d Sl --- f n 1 d Sn L p ( Sl )---L p ( Sn )QiLA ll A))(A 3l A))- 1 



o Jo 



(-1)" / ds! ■ ■ ■ ds„[X(s n )X(s n _i - s n ) ■ ■ ■ X(si - s 2 )X(t - si)]ij 
Jo 

-(-1)" / dfl!-.. /" \ Sn [Xs(0)Msn)Msn-l-S n )---X( Sl - S2 )±(t~ Sl )} lJ . (B5) 

JO Jo 
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By using this result, we can calculate 

^2(iLQe iLQs iLA lt A\)(Aj,A\)~ l = (X(t)) n - (X(0)X(t)) u 



- / d[(X(fl)X(t-a))ii-(X(0)X(s)X(t-«)) u ] 
Jo 

- / di / 1 ds 2 [(X( S2 )X( Sl - s 2 )±(t - si)) n - (X(0)X(s 2 )X( Sl - fl2 )X(i - si))n] 
Jo Jo 

+ • ■ • . (B6) 



In short, the Laplace transform of Eq. (|22|) is given by 



!Bi( l iIW TT^w- x<0)xL<s, rT^J n ' (B7) 



The other components of the matrix are calculated in the same way. 



[1] T. Koide, Phys. Rev. E75, 060103(R) (2007). 

[2] See, for example, P. Huovinen and P.V. Ruuskanen, Ann. Rev. Nucl. Part. Sci. 56, 163 (2006); Jean- Yves Ollitrault, Euro. 

J. Phys. 29, 275 (2008) and references therein. 
[3] I. Muller, Z. Phys. 198, 329 (1967); I. S. Liu, I. Muller and T. Ruggeri Ann. Phys. (N.Y.) 169, 191 (1986); as a review 

paper, see I. Muller, Living Rev. Relativ. 2, 1 (1999). 
[4] W. Israel, Ann. Phys. (N.Y.) 100, 310 (1976); W. Israel and J. M. Stewart, Proc. R. Soc. Lond. A 357,43 (1979); Ann. 

Phys. (N.Y.) 118, 341 (1979). 

[5] B. Carter, Proc. R. Soc. London, Ser A, 433, 45 (1991); as a review paper, see N. Andersson and G. L. Comer, Living 
Rev. Relativity 10, 1 (2007). 

[6] M. Grmela and H. C. Ottinger, Phys. Rev. E56, 6620 (1997); K. Dusling and D. Teanev. larXiv:0710.5932l 

[7] D. Jou, J. Casas- Vazquez, and G. Lebon, Rep. Prog. Phys. 51, 1105 (1988); 62, 1035 (1999). 

[8] T. Koide, G. S. Denicol, Ph. Mota, and T. Kodama, Phys. Rev. C75, 034909 (2007) . 

[9] G. S. Denicol, T. Koide, T. Kodama and Ph. Mota, Phys. Rev. C78, 034901 (2008). 
[10] G. S. Denicol, T. Koide, T. Kodama and Ph. Mota, J. Phys. G35 11 5102 (2008). 
[11] G. S. Denicol, T. Koide, T. Kodama and Ph. Mota. larXi"v:0808.3170l 

[12] W. A. Hiscock and L. Lindblom, Ann. Phys. (NY.), 151, 466 (1983); Phys. Rev. D31, 725 (1985); D35, 3723 (1987); 

Phys. Lett. A131, 509 (1988). 
[13] H. Kouno, M. Maruyama, F. Takagi and K. Saito, Phys. Rev. D41, 2903 (1990). 
[14] G. Torrieri and I Mishustin, larXiv:0805.0442l 

[15] T. Koide, G. Krein, and R. O. Ramos, Phys. Lett. B636, 96 (2006). 
[16] D. D. Joseph and L. Preziosi, Rev. Mod. Phys. 61, 41 (1989). 
[17] T. Koide and M. Maruyama, Nucl. Phys. A742, 95 (2004). 
[18] S. Nakajima, Prog. Theor. Phys. 20, 948 (1958). 
[19] R. Zwanzig, J. Chem. Phys. 33, 1338 (1960). 
[20] H. Mori, Prog. Theor. Phys. 33, 423 (1965). 

[21] As for review papers, see J. Rau and B. Muller, Phys. Rep. 272,1 (1996); U. Balucanai, M. H. Lee and V. Tognetti, Phys. 
Rep. 373, 409 (2003). 

[22] E. Fick and G. Sauermann, The Quantum Statistics of Dynamic Process (Springer- Verlag, Berlin, 1983). 
[23] L. E. Reichl, Modern Course in Statistical Physics, (University of Texas Press, 1980). 
[24] R. Zwanzig, Nonequilibnum Statistical Mechanics (Oxford Univerasity, New York, 2004). 
[25] T. Koide, Phys. Rev. E72, 026135 (2005). 

[26] R. Kubo, M. Toda and N. Hashitsume, Statistical Physics II (Springer- Verlag, Berlin 1983). 
[27] L. P. Kadanoff and P. C. Martin, Ann. Phys. 24, 419 (1963). 
[28] C. Cattaneo, Atti. Semin. Mat. Fis. Univ. Modena 3, 83 (1948). 

[29] P. M. Morse and H. Feshbach, Methods of Theoretical Physics (McGraw Hill, New York, 1953), p. 865. 

[30] K. Kawasaki, in C. Domb, M. S. Green (Eds.) Phase Transition and Critical Phenomena, Academic Press, London, 1976. 

[31] See, for example, D. R. Reichman and P. Charbonneau, J. Stat. Mech. P05013 (2005). 

[32] A. Yoshimori, Phys. Rev. E71, 031203 (2005). 

[33] K. Kawasaki, Physica A362, 249 (2006); J. Stati. Phys. 123, 711 (2006). 
[34] N. Hashitsume, F. Shibata, and M. Shingu, J. Stat. Phys. 17, 155 (1977). 
[35] F. Shibata, Y. Takahashi, and N. Hashitsume, ibid. 17, 171 (1977). 



15 



[36] F. Shibata and T. Arimitsu, J. Phys. Soc. Jpn 49, 891 (1980). 

[37] C. Uchiyama and F. Shibata, Phys. Rev. E60, 2636 (1999). 

[38] T. Koide, Prog. Theor. Phys. 101, 373 (1999). 

[39] T. Koide and M. Maruyama, Prog. Theor. Phys. 104, 575 (2000). 

[40] T. Koide, Prog. Theor. Phys. 107, 525 (2002). 

[41] T. Koide, Prog. Theor. Phys. 107, 1001 (2002). 

[42] J. Okada, I. Sawada and Y. Kuroda, J. Phys. Soc. Jpn 64, 4092 (1995). 

[43] This means that the memory function still has finite correlation after long time. This is a kind of long-time tail of memory 
function. This gives rise to a trouble when we interpret the equation as the Langevin equation because it is known that 
the memory function is given by the time correlation of the noise (The fluctuation-dissipation theorem of second kind). In 
this calculation, however, the noise term is just the component orthogonal to the gross variables and ignored. Thus, we 
will not discuss this problem in this paper. The physical meaning of the long-time tail may come from the incompleteness 
of the definition of the Mori projection operator. See, also the discussion in [17| ]. 

[44] Exactly speaking, there can exist the difference of the factor between the result of the POM and the GKN formula, because 
the linear response theory sometimes has an ambiguity for the definition of the transport coefficients. This is clear in the 
case of the definition of the diffusion coefficient. It is known that the POM can predict the correct expression, although 
the additional discussion is necessary to fix the proportional factor in the linear response theory [26| . 

[45] A. Hosoya, M. Sakagami and M. Takao, Ann. Phys. (N.Y.) 154, 229 (1984). 

[46] M. Ichiyanagi, Prog. Theor. Phys. 84, 810 (1990). 

[47] P. Kovtun, D. T. Son and A. O. Starinets, Phys. Rev. Lett. 94, 111601 (2005). 

[48] N. G. van Kampen, Physica 46, 315 (1970); W. L. Kath, Physica D12, 375 (1984); J. Dunkel, P. Talkner, and P. Hanggi, 

Phys. Rev. D75, 043001 (2007). 
[49] K. Kawasaki and J. D. G unton, Phys. Rev .A8, 2048 (1973). 
[50] H. Liu, D. Hou and J. Li, hep-ph/0609034 



